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Context. Most solar and stellar dynamo models use the a£l scenario where the magnetic field is 
generated by the interplay between differential rotation (the Q. effect) and a mean electromotive 
force due to helical turbulent convection flows (the a effect). There are, however, turbulent dy- 
namo mechnisms that may complement the a effect or may be an alternative to it. 
Aims. We investigate models of solar-type dynamos where the a effect is completely replaced by 
two other turbulent dynamo mechanisms, namely the Q. x J effect and the shear-current effect, 
which both result from an inhomogeneity of the mean magnetic field. 

Methods. We studied axisymmetric mean-field dynamo models containing differential rotation, 
the fl X J and shear-current effects, and a meridional circulation. The model calculations were 
carried out using the rotation profile of the Sun as obtained from helioseismic measurements and 
radial profiles of other quantities according to a standard model of the solar interior. 
Results. Without meridional flow, no satisfactory agreement of the models with the solar obser- 
vations can be obtained. With a sufficiently strong meridional circulation included, however, the 
main properties of the large-scale solar magnetic field, namely, its oscillatory behavior, its lati- 
tudinal drift towards the equator within each half cycle, and its dipolar parity with respect to the 
equatorial plane, are correctly reproduced. 

Conclusions. We have thereby constructed the first mean-field models of solar-type dynamos that 
do not use the a effect. 

Key words. Stars: magnetic fields - Sun: magnetic fields - magnetohydrodynamics (MHD) 

1. Introduction 

The standard dynamo model for the Sun and stars is the aQ model where, within the framework 
of mean-field magnetohydrodynamics, the magnetic field is produced by an interplay between dif- 
ferential rotation (the Q, effect) and the collective action of turbulent cyclonic convection flows, 
known as the a effect (Parker, 1955, 1979; Steenbeck et al., 1966; Krause & Radler, 1980). The a 
effect is here responsible for generating the poloidal component of the large-scale magnetic field 
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(LSMF), whose toroidal component is mainly generated by the the Q effect. The model is often 
supplemented with meridional flows, leading to so-called flux-transport dynamos (e.g., Choudhuri 
et al., 1995; Kiiker et al., 2001 ; Rempel, 2006; Dikpati & Oilman, 2007). The meridional flows may 
transport toroidal magnetic flux toward the equator and their speed may determine the cycle period, 
thus allowing us to bypass a number of problems connected with the a effect and aO, dynamos, 
as, for instance, that, in the case of the Sun, the obtained cycle periods are generally too short and 
the magnetic activity is not sufficiently concentrated at low latitudes (see, e.g., Ossendrijver, 2003; 
Riidiger & Hollerbach, 2004; Brandenburg & Subramanian, 2005a). 

In mean-field magnetohydrodynamics, the influence of the turbulence on the LSMF is ex- 
pressed by the mean turbulent electromotive force (MEMF), £ = (m x b), where u and b are 
the fluctuating parts of the velocity and magnetic field and angular brackets denote averages. The 
by far best known contribution to £ is provided by the a effect, namely, a turbulent electromotive 
force a{B}, with a denoting a (symmetric) tensorial factor of proportionality and {B} the LSMF. 
However, there are other turbulent dynamo mechanisms besides the a effect. Two of them are the 
Q X / effect (Radler, 1969; Stix, 1976) and the shear-current or W x / effect (Rogachevskii & 
Kleeorin, 2003, 2004); SI is here the angular velocity of the stellar rotation, J - W x {B)I^q the 
large-scale electric-current density, and W - V xV the large-scale vorticity, V denoting the large- 
scale velocity. Both these effects result from an inhomogeneity of the LSMF, in contrast to the a 
effect, which also works with a homogeneous {B) (that is to say, for calculating the a effect, {B) 
may be considered as homogeneous on the scale of the fluctuations). 

In the commonly used representation of the MEMF on the basis of symmetry considerations 
(see Radler, 1980; Krause & Radler, 1980; Radler, 2000; Radler et al., 2003), the Q x / and shear- 
current effects represent contributions to the S term, a term of the form £,5 = x (V x {B)), where 
5 is a vector. Since £>s J = 0, the effects described by this term cannot bring energy into the 
mean magnetic field and, thus, cannot lead to working dynamos when acting alone. These effects 
have been investigated little in the context of solar and stellar dynamos. For a recent study of the 
possible role of the fix J effect when acting together with the a effect and differential rotation in 
a spherical shell, or when acting together with another part of the MEMF, not included in dynamo 
studies before, in a rigidly rotation full sphere, see Pipin & Seehafer (2009), where an illustration 
of the physical mechanism behind the il x J effect also may be found; the mechanism of the 
shear-current effect is very similar to that of the Six J effect. 

In this paper, we consider mean-field dynamo models in the geometry of a spherical shell, 
as appropriate for solar-type stars, where the a effect is completely omitted. Instead, the Six J 
and shear-current effects serve as turbulent dynamo mechanisms. In nearly all mean-field dynamo 
studies, the effective strengths of the different physical ingredients are controlled by freely varied 
dimensionless parameters; in the case of dynamo effects, e.g., the a effect, these are usually re- 
ferred to as dynamo numbers. This reflects our present knowledge of the physical processes in the 
convection zones of the Sun and stars. Realistic self-consistent numerical models of these processes 
and their interactions will remain out of reach for the foreseeable future. Given this situation, we 
deem it advisable to explore the potentials of turbulent dynamo effects other than the a effect. 

Numerical evidence for turbulent dynamo effects has so far mainly been obtained from con- 
vection simulations in small (compared to the dimensions of a star) rectangular boxes (e.g., 
Brandenburg et al., 1990; Ossendrijver et al., 2001, 2002; Giesecke et al., 2005; Kapyla et al.. 
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2006; Cattaneo & Hughes, 2006; Hughes & Cattaneo, 2008). Due to the assumption of a uniform 
mean magnetic field and other hmitations, most of these studies could only find parts of the MEMF 
that are proportional to the LSMF, i.e., the a effect and turbulent pumping (a contribution to the 
MEMF of the form £y = 7 x (B), with 7 denoting a vector; it leads to an advection of the mean 
magnetic field). Recently, however, Kapyla et al. (2009), who used a procedure referred to as the 
test field method (Schrinner et al., 2005, 2007) together with numerical simulations of turbulent 
convection with shear and rotation, were able to also identify the action of the combined Six J and 
shear-current effects. 

Here, we explore axisymmetric kinematic dynamo models containing the SI x J and shear- 
current effects, differential rotation, and meridional circulation. In calculating the MEMF we use 
analytical expressions derived by Pipin (2008) on the basis of a simplified version of the r approx- 
imation (cf. Vainshtein & Kichatinov, 1983; Brandenburg & Subramanian, 2005a,b). We construct 
models with distributed dynamo action in the bulk of the convection zone, rather than in the over- 
shoot layer at the bottom of the convection zone. The model calculations are carried out using the 
rotation profile of the Sun as obtained from helioseismic measurements and radial profiles of other 
quantities according to a standard model of the solar interior 

The remainder of the paper is organized as follows: In Sect. 2 we describe our dynamo model, 
as well as the used numerical procedure (some benchmark tests for our computer code are presented 
in Appendix A). Then, in Sect. 3, we present the obtained results. In Sect. 4, we draw conclusions 
and discuss our results. 



2. Model and numerical procedure 

The axisymmetric LSMF is written in the usual way as the sum of a poloidal and a toroidal part, 

{B}^Vx—^ + Be^, (1) 

where A(r, 9, t) (the flux function for the poloidal field) and B{r, 0, t) (the toroidal field component) 
are scalar functions of radius, r, colatitude, 0, and time, t, and is the unit vector in the direction 
of the azimuthal coordinate, (p. The mean-field induction equation then takes the form 
dA . UgdA OA 

— ^rsm96^ - U,-— , (2) 

ot r 06 or 

dB \d(Sl,A) ^\l dr(&e-UrB) 8 (5, + UgB) \ 



dt r d{r,e) r\ dr 86 
where the effects of the large-scale flows enter through the differential-rotation rate, Q(r, 9) = 
\Sl{r, 0)\, and the components of the meridonal flow, {/, and Ug. 

To calculate the MEMF, whose effects appear through the components of £ in Eqs. (2) and (3), 
we modify the expressions given in Pipin & Seehafer (2009) by completely omitting the a effect 
but including, in addition to the Six J effect, isotropic and anisotropic turbulent diffusion, and 
turbulent pumping, now also the shear-current effect. The contribution of this to the MEMF is to 
linear order, i.e., for a weak mean magnetic field, as well as neglecting the effect of the Coriolis 
force, given by (Pipin, 2008) 

= Sinm\C\VlmB„i + ClBnlVml + CiVlmBin 

9\ (4) 
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where tensor notation and the summation convention have been used. Bij - d{Bi}/dxj is the 
gradient tensor of the mean magnetic field and Vij a corresponding quantity for the differential 
rotation, namely, V/y = dVi/dxj, with V = r sin (f2 - Qo) being the rotational velocity, de- 
fined in a reference frame that rotates with angular velocity Qq, the rigid-body rotation rate of the 
core (encountered at midlatitudes through the convection zone). is a small-scale or turbulent 
convective background velocity as present in the absence of rotation and a mean magnetic field, 
Tc the correlation time of m*°', and Ci = (e- 3/5) /6, C2 = (e- 1) /5, C3 = (1 -H e) /15, and 
C4 - - (7e +11) /30 are constants; s — ^m^, y/JIopj is the square root of the prescribed 

ratio between the energies of a fluctuating magnetic background field 6*"', assumed to be generated 
by a small-scale dynamo (which is fully independent of the mean magnetic field), and the back- 
ground velocity field h'°* {u^ - -^^h*"'^^ is the rms value of the convective background velocity 
field and p the mass density). 

In the foUowing, we assume energy equipartition between the two background fields, i.e., e = 1 . 
Furthermore, only the azimuthal Six J and shear-current effects are taken into account. This may be 
justified by the fact that the toroidal part of the solar LSMF is much stronger than the poloidal one. 
However, the remaining parts of the MEMF (isotropic and anisotropic turbulent diffiision, turbulent 
pumping) are included in all components. The components of the MEMF in spherical coordinates 
then become 



{ ff + Iff" sin^ 6 d sin 9B sin 26 drB 

1 rsind 86 r dr 



+ G sin 26 ff^B 



ff + 2ff cos^ 6 drB 2ff cos 6 8 sin 6B 
r dr r do 

-G(/f^+cos2<')B|, 



(5) 



ff+2f["^ I I d'^A I d I dA 

&^^f,T\ : \—n^ + 



(6) 



G 



sm6 dr^ r^-dOsindde 
2 cos 6 &^dA ff" cos^ 6 +&> dA] 



■2 do rsinO dr 

<n) Ad) 



(7) 



■'^ \ dr r d6 

(p ' 

with £^*^' in Eq. (7) denoting the contribution of the shear-current eff'ect, given by 

11 /A .\/ dB sine dB 



d6 



1 5 b,q) 

H- -sine^:^ 

3 d(r,e) 

Bl dCl sin6dQ. 
+ — [cos 6 



(8) 



6 \ dr r do 

ff \ ff\ ff \ and /J''' denote functions of e and the Coriolis number Q* = 2Q.oTc that are given 
in Appendix B, G = {d/dr)logp is the density scale factor, Q - Q./Q.0, and fir = Cr^rjT, with 
C,, < 1, C^"' < 1, C^*^' < 1 are parameters to control the relative strengths of 
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different turbulence effects. C,, regulates the turbulence level, and C^"' and C^'*^' weight the ilx J 
effect and the shear-current effect, respectively. 

Currently, the dependence of the shear-current effect on the Coriolis number is unknown. 
Eq. (4) has been derived disregarding the effect of the Coriolis force and is, thus, safely appli- 
cable only in the limit of slow rotation, Q* «c 1. But in the solar convection zone, in particular its 
deeper layers, the Coriolis number is large, Q* » 1 (cf., e.g.. Fig 2 in Pipin & Seehafer, 2009). To 
take this into account, we modulate the value of Sj)*^*, given by Eq. (8), by the quenching function 
f^''^ (Q*) which also appears in the expression for the Six J effect (penultimate term in Eq. (7), 
proportional to Cf^^); additionally, the expression for the shear-current effect is normalized so as 
to give Eq. (4) in the limit of slow rotation, for which one finds /^'" (Q*) ~ (1/5) n*. Without such 
a quenching, i.e., directly applying Eq. (4), the shear-current effect would become unrealistically 
strong at the bottom of the convection zone. 

In our numerical calculations we have used a dimensionless form of the equations, substituting 

r - xRq and t rjot/R^, where rjo - 1.8 ■ 10^ m^/s is the maximum value of tjt in the convection 

zone; that is, length is measured in units of the solar radius and time is measured in units of the 

turbulent magnetic diffusion time. To ~ 8.5 yr. The integration domain is radially bounded by 

X = 0.72 and x - 0.96. The boundary conditions on the magnetic field are the usual approximate 

dxB 

perfect-conductor conditions, i.e., A — 0, — 0, at the bottom boundary (Kohler, 1973), and 

OA- 
vacuum conditions, that is, B = and a continuous match of the poloidal field component to an 

exterior potential field, at the top boundary. 

The radial profiles of characteristic quantities of the turbulence, such as the rms value and 

the coiTelation length and time, £c and Tc, of the convective background velocity field, as well as 

the density stratification, were calculated on the basis of a standard model of the solar interior (Stix, 

2002), assuming the ratio of the correlation length to the pressure scale height (referred to as the 

mixing-length parameter) to be 1.6. The rotation profile as known from heUoseismic measurements 

(Schou et al., 1998) is approximated by 

Q = Qo/(a,0), (9) 

with 

fix, 9) =^ [435 +51 {x - xo) + 264>(x) (l - 5 cos' o) (10) 
-3.5(1 - 14cos2 6> + 21cos''6l)] , (11) 

where 

(^(jc) = 0.5{1 +tanh[50(ji:-xo)]} (12) 

and X() = 0.71 is the position of tachocline, situated below the bottom boundary of the integration 
domain. 

A remark concerning our locating the lower boundary at x = 0.72 seems in order. Very often 
this boundary is placed at x = 0.65 (cf., e.g., Jouve et al., 2008). Then, however, some modeling 
of the tachocline is needed, where the differential rotation changes into rigid rotation in the ra- 
diative core. The physical parameters of this transition region are rather uncertain at the moment. 
Here, we consider a convection-zone dynamo model with distributed dynamo action in the bulk 
of the convection zone, where all physical parameters needed can be derived from helioseismic 
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measurements and the standard model of the solar interior. There are other dynamo models where 
the dynamo just operates in the tachocline (a critical discussion of arguments for and against deep- 
seated and distributed dynamos is found in Brandenburg, 2005). 

The meridional flow, U, is modeled in the form of two stationary circulation cells, one in the 
northern and one in the southern hemisphere, with poleward motion in the upper and equatorward 
motion in the lower part of the convection zone (for the theory of the meridional circulation see, 
e.g., Rempel, 2005, 2006; Miesch et al., 2008; Brun & Rempel, 2009). The condition of mass 
conservation, V ■ (pU) - 0, is ensured by a stream-function representation of pU (cf., e.g., Bonanno 
et al., 2002), so that 

t7=-Vx-;-^ = -V^x— ^. (13) 
p xsmO p xsmO 

The stream function, tfr, is written as 

(A(x,6/) = MoiAfcf), (14) 

with Mo denoting the maximum amplitude of U, which is treated as a free paramter. Our choice for 
(^(x, 6) is 



4/{x, 0) =co sin^ cos T(x) , 



CO 



max 



1 sin 0cos,6T{x)ej 
V X 



(15) 



p jc sin 

where the function T{x) is selected such that the top and bottom boundaries are impenetrable and 
stress-free (free of tangential stresses), that is (see, e.g., Batchelor, 1967), 

t^x = 0, -t- - [/fl = atx = xt, x = Xb, (16) 

o0 ox 

where Xt and Xb are the outer and inner radius, respectively, of the considered spherical shell. Using 

Mo Co (3 cos^ 0-\)t uq Co sin 6* cos dT 

U,^ ^ 2 Ue^-— — , (17) 

px^ px dx 

as follows from Eqs. (13), (14), and (15), the boundary conditions given by Eq. (16) take the form 

d 1 dr 1 dr 

r = , X =0 al X - Xt, X — Xb (18) 

dx px dx px dx 

and are satisfied with 

r (X) = P2 (^) + Ci [P3 (f) - ^] + C2 [P4 (f) - 1] - 1 , (19) 

where 

2x — (xt + Xh) 

^ = — ^ (20) 

— Xb 

is the radial coordinate transformed from the interval [xb, Xt] to the interval [-1, 1], P„ denotes the 
Legendre polynomial of degree n, and ci ^ -0.207 and C2 ~ -0.097 are numerically determined 
constants; for the normalization constant cq (cf. Eq. (15)) one then finds cq ~ 0.027. 

Fig. 1 shows radial profiles of the mass density, the isotropic and anisotropic magnetic diffii- 
sivities, and the effective strength of the £lx J effect, and in Fig. 2 the effective strengths of the 
contributions to the shear-current effect as given by Eq. (8) are displayed. The properties of the 
two large-scale flows used in the model, namely, contours of the differential rotation rate and the 
geometry and strength of the meridional circulation, are depicted in Fig. 3. Other model quantities, 
such as the radial profile of the Coriolis number, Q*, and the pumping velocity of the toroidal part 
of the LSMF, are as in Pipin & Seehafer (2009, see Fig. 2 there). 
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Fig. 2. Effective strengths of the contributions to the shear-cuiTent effect as given by Eq. (8): (11 /6) (Cl- l^N 
(top left), {l/6){cosedCl/dx-(smd/x)dCl/dejN (top right), (l / 3) sin d {dCl / dOj N (bottom left), and 
-(1/3) sine [dCl/dr) N (bottom right), with N = /f' (ft*) / (/f' (O*) + ifi"^ (fi*)). 



Numerical procedure 

In our numerics we use a Galerkin method, expanding the magnetic field in terms of a basis that 
satisfies the boundary conditions implicitly. The system of Eqs. (2) and (3) admits exponentially 
growing or decaying solutions, which we represent in the form 

A (x, e, t) = e-^^ ^ ^ A„„, sin e^*-^,' (^) P]„ (cos 6) , (21) 

n m 

B {X, 0, t) ^^^^YjYj ^'""^ ' ^22^ 

n m 

where S^^,l and are linear combinations of Legendre polynomials, and f is the associated 

Legendre function of degree m and order 1 . These expansions ensure the regularity of the solutions 

at the poles 6 - Q and 9 - n, where Bq and Ba,, i.e., and B, have to vanish. The radial boundary 

sinfl 

conditions are satisfied by the choice (see "basis recombination" in Boyd, 2001; Livermore & 
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Fig. 3. The large-scale flows. Left panel: Contours of the rotation rate Q (left) and streamlines of the meridional 
flow, U, i.e., contours of its stream function, iff ( right). Right panel: Radial profile of the component Ue of the 
meridional flow at a latitude of 45° (in m/s). 



Jackson, 2005) 



S^tJl (^) = P„-i + a^Pn (^) + flr«;n+i (^) , 



where 



2n + l 



(n + 1)2 + 2m/y ' 



xb(« + 1)2 - 2/r 



n + 2m I y 
(n + 1)2 + 2m/y ' 
Xb»2 - 2/y 

xb(«+ i)2-2/r ' 



(23) 
(24) 

(25) 
(26) 



with y - 



denoting the derivative of ^ with respect to x. 



Integrations over radius and latitude, necessary for calculating the expansion coefficients a„,„ 
and b„,„, were done by means of the Gauss-Legendre procedure, and the eigenvalue problem for 
determining the exponent A and the associated eigenmodes was solved by means of Lapack rou- 
tines. The spectral resolution was 15 modes in the radial basis and 22 modes in the latitudinal 
basis for the calculations of growth rates (including stability boundaries) and dynamo periods, and 
14 X 30 modes for simulations of time evolutions and butterfly diagrams (see Sect. 3 below); by the 
assumption of either dipole-type or quadrupole-type solutions the latitudinal resolution could be 
doubled in a part of the calculations. The results were qualitatively confirmed by a number of runs 
with still higher resolution. Benchmark calculations for the code used are presented in Appendix 
A. 



3. Results 

3.1. 5*"'Q dynamo with meridional flow 

The meridional flow becomes essential for the dynamo if the effective magnetic Reynolds num- 
ber, based on the meridional flow velocity and the turbulent magnetic diffusivity, is high enough. 
Meridional flow velocities higher than about 10 ... 20 m/s can scarcely be brought into agreement 
with the solar observations. Thus, the turbulent magnetic diffusivity should be low. In our formu- 
lation, all turbulence effects are consistently scaled by the parameter Crj. Decreasing Crj leads to 
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Fig. 4. Difference between the growth rates of the first (most unstable) dipolar mode and the first quadrupolar 
mode in the plane spanned by C^^^ and uq for C^*^' = 0. Red colors correspond to a dominance of dipolar 
modes and blue colors to a dominance of quadrupolar modes. Also shown are isolines of the frequency of the 
first dipolar mode (in areas where it oscillates), with numbers giving the oscillation frequency {In IT, T being 
the time period of the oscillation) in units of the inverse magnetic diffusion time, t/q/R^. The solid bold line 
indicates the stability boundary for the first dipolar mode, with the unstable (dynamo) region lying to the right 
of this line. 



increasing the influence of the flow on the magnetic field and acts, thus, like increasing the am- 
plitude of the flow. An advection-dominated regime with a solar-like magnitude of the meridional 
flow is obtained if C;, g 1 /20. Below, we fix C,, to the value 1 /40. 

Figure 4 illustrates the bifurcation scenario for a dynamo with meridional flow. Here 

the sources of the dynamo are the £lx J effect (the shear-current effect being neglected) and dif- 
ferential rotation. Also included are turbulent diffusion (isotropic and anisotropic) and turbulent 
pumping. As the oberservations of solar activity suggest, the large-scale solar magnetic field is 
characterized by an antisymmetric parity with respect to the equatorial plane. In the figure, the 
plane spanned by the parameters C^"' (measuring the strength of the ilxj effect) and uq (the max- 
imum amplitude of the meridional flow) is displayed, with red/blue colors indicating a dominance 
of dipolar/quadrupolar modes, which are antisymmetric/symmetric with respect to the equatorial 
plane, over quadrupolar/dipolar modes. For a sufficiently high velocity of the meridional flow the 
dipolar parity dominates, as needed for the Sun. The critical dynamo number cf^^ where the first 
dipolar mode becomes unstable does not depend very much on the flow strength. This may result 
from the fact that the meridional flow mainly acts as a conveyor belt for the magnetic field, rather 
than as a generation mechanism. 

In Fig. 5 (solid line) the dependence of the dynamo period on the amplitude of the meridional 
flow along the stability boundary is shown. As expected, the period is a decreasing function of the 
flow velocity (cf. Bonanno et al., 2002, 2006). The dependence approximately follows a power law 
with a scaling exponent of -0.7. 

Figure 6 (top and middle) shows, in the form of the strength of the toroidal field and field lines 
of the poloidal field in the meridional plane, the evolution of the LSMF in an dynamo model 
on the basis of the first dipolar mode for C^"^ — 0.006 and uq - 15m/s. In the bottom panel 
of the figure, an associated simulated butterfly diagram, i.e., the strength of the toroidal LSMF 
(integrated over depth in the convection zone) in the time-latitude plane is shown, together with 
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Uq, [m/s] 

Fig. 5. Dependence of the dynamo period on uq along the stability boundary of the most unstable dipolar mode 
for the (5*"'n (solid line) and S'-^^Q. (dashed line) dynamos. 

isocontours of the radial field component at the top boundary of the considered spherical sheU. 
The latitudinal drift of the toroidal field towards the equator, as indicated by the observations of 
sunspots, is qualitatively correctly reproduced; the poloidal field is coupled to the toroidal field 
and shows a similar drift towards the equator. The phase relation between the toroidal and poloidal 
parts of the magnetic field, in particular the polar reversal of Br shortly after the maximum of 
at low latitudes, with BrB^ < before and BrB^ > after the polar field reversal, is also in good 
agreement with the solar observations. 

In the example shown in Fig. 6, the obtained cycle period is about 80 yr, which is nearly 
four times the period of the solar activity cycle. Tuning the paramters cannot significantly reduce 
the period. This may appear not fully satisfactory, but the period is at least in the right order of 
magnitude. Furthermore, for the example in Fig. 6, we find B^/B^ ~ 15, which is by a factor of 
10 ... 100 smaller than the ratio between the large-scale toroidal and poloidal fields as supposed 
for the solar convection zone. In general, increasing the speed of the meridional flow in the model 
reduces the obtained ratio between the toroidal and poloidal fields, in apparent conflict with the 
need to reduce the cycle period. 

3.2. dynamo with meridional flow 

Figures 7 and 8 show results for a dynamo on the basis of the shear-current effect, differential ro- 
tation, and a meridional flow, with the a effect and the £lxj effect being neglected (but again with 
turbulent diffusion and turbulent pumping being included). The results strongly resemble those for 
the 5*"*f2 dynamo with meridional flow as described in Sect. 3.1, obviously due to the similarity 
between the Q x / effect and the shear-current effect. The comments given in Sect. 3.1 thus ap- 
ply here as well. The dependence of the dynamo period on the amplitude of the meridional flow 
along the stability boundary (Fig. 5, dashed line) approximately follows a power law with a scaling 
exponent of -0.6. 

4. Conclusions 

We have studied kinematic axisymmetric mean-field dynamo models in the geometry of a spherical 
shell, as appropriate for the Sun and solar-type stars, where the il x J and shear-current effects 
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TIME, [Yr] 

Fig. 6. S^^'Q dynamo for C^"' = 0.006 and uq = 15m/s. Top and middle: Snapshots of the strength of the 
toroidal LSMF (color-coded, red/blue colors correspond to positive/negative field values) and field lines of the 
poloidal LSMF (solid/dashed lines indicate clockwise/counter-clockwise field direction) over half a cycle in 
the meridional plane. Bottom: Associated butterfly diagram, showing the (color-coded) strength of the toroidal 
LSMF (integrated over depth in the convection zone) in the time-latitude plane. Also shown are isocontours of 
the radial field component at the top boundary (solid red/dashed blue lines indicate positive/negative values). 




Fig. 7. As Fig. 4, but with replaced by Cj", and = 0. 

were included as turbulent sources of the large-scale magnetic field while the a effect was omitted. 
Besides the turbulent dynamo mechanisms and differential rotation, a meridional circulation, in the 
form of two stationary circulation cells, one in the northern and one in the southern hemisphere, 
also was incorporated into the models. We have concentrated on the dynamo onset. 
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50 100 150 

TIME, [Yr] 

Fig. 8. As Fig. 6, but for a (J^^'fi dynamo with C^' = 0.007 and uq = 15 m/s. 

Our results show that the Six J and shear-current elfects can, at least in principle, take over the 
role that the a effect usually plays in mean-field dynamo models. However, only if the meridional 
flow is sufficiently fast are the characteristic properties of solar-type dynamos qualitatively cor- 
rectly reproduced. In particular, the amplitude of the meridional flow needs to exceed a threshold 
value in order that the most unstable magnetic mode has dipolar parity and oscillates. This mode 
then also shows a latitudinal drift towards the equator within each half cycle and a phase relation 
between the poloidal and toroidal parts of the field in accordance with the observations of solar 
activity. The threshold value for the amplitude of the meridional flow, uq > 10 m/s (reached at the 
surface, the flow speed at the bottom is on the order of 1 m/s), is consistent with solar observations 
and agrees with the value of 10 m/s often adopted in studies of flux-transport dynamos with the a 
effect (cf., e.g., Bonanno et al., 2002; Jouve et al., 2008). 

In models of advection-dominated dynamos, the specifics of the turbulent dynamo mechanism 
that generates the mean poloidal field are less important than they are in models without merid- 
ional flow (given the rotation law and the generation of the mean toroidal field from the mean 
poloidal field by velocity shear). Once the field is generated, it is advected equatorwards by the 
flow. However, the distribution of the turbulent dynamo sources, or their more or less strong local- 
ization, decisively influences the parity properties of the LSMF. Studies of flux-transport dynamos 
with an a effect as the turbulent source of the LSMF indicate that the a effect must be strongly 
localized at the bottom of the convection zone to ensure the correct (dipolar) parity of the LSMF 
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(Dikpati & Gilman, 2001; Bonanno et al., 2002). In our models, the turbulent dynamo effects are 
distributed over the bulk of the convection zone, though they are strongest near the bottom of the 
included domain. We note that our turbulent dynamo sources are not introduced arbitrarily but are 
calculated using a standard model of the solar interior together with rotation rates obtained from 
helioseismic measurements. 

As other advection-dominated dynamo models, the models presented here work only if the 
effective magnetic difFusivity is strongly reduced compared to the mixing-length estimates. At a 
radial distance of, say, 0.85 Rq, the turbulent magnetic difFusivity in our models is about 0.2 770 {t]o - 
1.8 ■ lO'^m^/s is the maximum value of the turbulent magnetic difFusivity in the convection zone 
according to the mixing-length estimate). Together with our value of 0.025 for the parameter Crj 
(which regulates the turbulence level), this gives an effective magnetic diffusivity of about 10^ m^/s, 
in agreement with the upper limit of 3 ■ 10^ m-/s for the turbulent magnetic diffusivity in the bulk 
of the convection zone given by Dikpati & Gilman (2006, 2007) for flux-transport dynamos with 
the a effect. 

The cycle periods that we obtain are at least three times as long as the observed period of 
the solar activity cycle. Also, the ratio between the toroidal and poloidal parts of the large-scale 
magnetic field is significantly smaller than supposed for the solar convection zone (apparently a 
common problem of all flux-transport models). Here one should keep in mind that requiring a 
perfect fit to the solar details, as far as these are known, would overstress the models. For instance, 
a solution that bifurcates at the dynamo onset will change quantitatively if it is traced away from 
the bifurcation point. Thus, ultimately, self-consistent nonlinear models will be needed. 
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Appendix A: Benchmarks for the code used 

Here we present some benchmark tests for the computer code that we used and show how the 
crititical dynamo numbers and dynamo periods for the models of Sects. 3.1 and 3.2 converge for 
an increasing number of modes taken into account. 



A. 1 . Free decay modes 

First we test the accuracy of the implementation of the exterior boundary conditions and the speed 
of convergence. If we neglect all the dynamo effects in Eqs. (2) and (3), only simple isotropic 
diffusion remains and the equations take the form 

dA _ d^A sine d 1 dA 

'dt ~ ^~'d9'^~de' '^^'^^ 

dB _l d^ jxB) 1 d 1 g (sin eg) 
'dt ~ ~x~d^ l^-'de'^ m ' 

where for simplicity the magnetic diffusivity has been assumed to be homogeneous and has been set 
equal to unity. The equations for the poloidal and toroidal parts of the field are decoupled here. The 
eigenmodes to Eqs. (A.l) and (A. 2) are the free decay modes, exponentially decaying oc exp/l,f, 
where the /I,- are the eigenvalues of the Laplacian operator for the considered domain under the 
imposed boundary conditions; these eigenvalues are all real and negative. For the test, we consider 
the case of a full sphere (rather than a spherical shell) surrounded by vacuum, for which the free 
decay modes can be determined analytically and are well documented in the literature (see, e.g., 
Moffatt, 1978; Backus et al., 1996). For that purpose, the potential functions A and B are written as 

A (^, e, t) ^^^'YjYj ^«"' ^ ^ ™ ' ^^-^^ 



n m 



B (x, 6, t) = e-^' 2 2 B„„S ix) pI (cos 0) , (A.4) 



with 



cWr^ /p n r , (2« + l)(2« + 2) + 2(m+l) \ 
S„m W = I ^2,7+1 (x) - Pi (x) 2mT4 1 ' 

5^^) (x) = X (P2„+i (x) - Pi (x)) , (A.6) 

where the radial variable, x, varies in the interval [0, 1]; the transformation to the variable ^ (cf. Eqs. 
(21) and (22) in Sect. 2) is not used here. By the choice of the basis functions given by Eqs. (A. 5) 
and (A.6) the exterior vacuum conditions are satisfied and the regularity of the fields at the origin 
is ensured. This set of basis functions differs from that used in our calculations for the spherical 
shell, but the general structure of the code and the solution algorithms are not changed. 

The dependence of the solutions to Eqs. (A.l) and (A. 2) on radius is given analytically in terms 
of the spherical Bessel functions j„{x) oc ^1/ V^) y„+i/2(x) (where 7„+i/2 is the ordinary Bessel 
function of half-integer order n + I /2), and the decay rates, -A,,, are given by the squares of the 
zeros of the functions y„_i for the poloidal and y„ for the toroidal modes. The slowest decaying 
(largest scale) poloidal mode decays with the rate = n-, the slowest decaying toroidal mode 
with the rate ^ 4.4934^; the corresponding eigenfunctions are S'^^\x,ff) oc ji (^a'-^^ x^ sin^ 

and S ^f\x, 0) oc ji (^Af^x^ sin 0. These two decay modes are used for the test. The dependences of 
their potential functions are given by the first terms (with m - 1 ) in the latitudinal expansions on the 
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Convergence of the eigenvalues and eigenvectors of the slowest decaying poloidal and toroidal 



N 


E(A) [A] 


E{B) [A] 


E(A) [B] 


E{B) [B] 


3 


1.08e-7 


3.98e-9 


3.83e-5 


6.849e-4 


4 


5.651e-ll 


1.255e-12 


9.984e-8 


4.365e-9 


5 


9.68e-14 


2.142e-16 


1.207e-10 


2.497e-12 


6 


7.72e-14 


2.136e-20 


8.707e-14 


9.068e-16 


7 


1.38e-14 


1.325e-24 


1.77e-14 


4.80e-19 


8 


1.90e-14 


4.427e-29 


5.32e-15 


6.263e-23 



^ A' is the number of modes in the radial basis, and E(A) [A] and E(B) [A] are the errors of eigenvalue and 
eigenvector for the poloidal mode and E(A) [B] and E(B) [B] the corresponding errors for the toroidal 
mode. 

right-hand sides of Eqs. (A. 3) and (A.4). (The potential functions A and B differ from the potentials 
S and T in the poloidal-toroidal decomposition B = V x (r x WS) + r x WT as normally used in 
non-axisymmetric cases, see, e.g., Moffatt (1978) and Backus et al. (1996). For our axisymmetric 
case, one has dS jdO - A/ sinO and dT/dO - B. The angular dependence of both S for the slowest 
decaying poloidal mode and T for the slowest deacaying toroidal mode is given by the spherical 
surface harmonic ^^(cos 0) oc cos 9, in agreement with the 6 dependences of A and B as given 
above.) 

Table A.l shows the convergence of the eigenvalues and of the corresponding eigenvectors 
for our numerical scheme. Similar to Livermore & Jackson (2005), the eigenvectors are scaled 
so that 5 j'*' (x - 1,0 - 7t/2) - 1 and 5 j^' (x - 0.5, — 7r/2) — 1, and the errors are measured as 
E (A) — |/iti-ue - Anuml and E (B) — Jy |Btrue " ^numl" dV. The number of modes in the radial basis, 
A^, is varied, while in the latitudinal basis just the first mode is taken into account. The convergence 
is seen to be exponential in both the poloidal and toroidal cases. 

A.2. Test case B ofJouve et al. (2008) 

The next test case is taken from Jouve et al. (2008), who presented a comparitative benchmark 
study of different numerical codes for axisymmetric mean-field solar dynamo models in spherical 
geometry. Here we consider their test case B, which is a pure aQ. dynamo in a spherical shell with 
sharp gradients of the turbulent magnetic diffusivity and the strength of the a effect at the bottom 
of the convection zone; for details we refer to Jouve et al. (2008). The potential functions A and 
B are expanded according to Eqs. (21)-(26) in Sect. 2, and the integration domain is now radially 
bounded by x^ = 0.65 at the bottom and Xt = 1 at the top. 

In Jouve et al. (2008), the strength of the a effect is regulated by a dynamo number, €„■ The 
different codes are compared by indicating in tables the critical a-effect dynamo number, C™', at 
which exponentially growing solutions appear, and the corresponding oscillation frequency, a> - 
2n/T. In addition, butterfly diagrams and the evolution of the fields in the meridional plane are 
shown. Our values of C™' and oj for different spectral resolutions are given in Table A. 2, and 
Fig. A.l shows the temporal evolution of the toroidal and poloidal parts of the field (i.e., of the 
unstable eigenmode) at the critical dynamo number The values in Table A. 2 are in best agreement 
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Table A.l. 

modes. " 
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Table A.2. Test case B of Jouve et al. (2008) ^ 
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Resolution C"'' o) 



8x8 


0.443 


180.5 


10 X 10 


0.4175 


175.1 


12 X 12 


0.4095 


172.2 


13 X 13 


0.411 


172.4 


14 X 14 


0.4122 


172.7 


16 X 16 


0.4125 


172.9 



^ Critical values of the dynamo number, C"", and the oscillation frequency at the dynamo onset, ai, for 
different radial and latitudinal spectral resolutions, co is measured in units of the inverse magnetic diffusion 
time. 







\=5-nJ6 



Fig. A.l. As Fig. 6 (top and middle), but for test case B of Jouve et al. (2008) at the critical a-effect dynamo 
number with a spectral resolution of 16 x 16 modes. 



with those given in the corresponding table, Table 3, of Jouve et al. (2008). Similarly, the evolution 
shown in Fig. A.l is apparently identical to that shown in the corresponding figure. Fig. 7, of Jouve 
et al. (2008); the same applies to the simulated butterfly diagrams (not shown here). 

A.3. Convergence of crititical dynamo numbers and dynamo periods for the models of 
Sects. 3. 1 and 3.2 

Fig. A. 2 shows the convergence of the critical dynamo numbers (where the first dipolar mode 
becomes unstable) and associated dynamo periods for the 5*"'f2 dynamo model considered in 
Sect. 3.1 and for the 5**^'f2 dynamo model considered in Sect. 3.2. The amphtude of the meridional 
flow is Mo - 15 m/s and uq - 25 m/s; uq - 15 m/s is the value we used most, and uq - 25 m/s is 
the highest meridional-flow amplitude that we considered, corresponding to the largest magnetic 
Reynolds number in the study. High Reynolds numbers are known to cause numerical problems. 
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Fig. A.2. Convergence of the critical dynamo numbers (left) and associated dynamo periods (right) for the 
models of Sects. 3.1 and 3.2. A' is the total number of modes taken into account. Calculations were done for 
resolutions of 6 x 6, 6 x 8, 8 x 10, 10 x 14, 16 x 18, 15 x 22, and 14 x 30 modes in the radial and latitudinal 
bases, respectively (in addition, the dipolar symmetry was taken into account, so that the highest latitudinal 
resolution is actually 50). Solid lines refer to the i5'"*n dynamo model and dashed lines to the i5"''*n dynamo 
model, and blue color corresponds to iig = 15 m/s and red color to ug = 25 m/s. 



Appendix B: Definitions of the functions f."^ and /]' 

Here we give the definitions of the functions f^"^ and f^'^^ that are used in the representation of the 
turbulent electromotive force £. For details of the calculations we refer to Pipin (2008). 

1 



A") _ 

■13 



■12 



Ad) 
■I A 



4n*2 

1 

4^*2 
1 

4Q*2 

-(3e+ 1) 
1 



/ . \ arctan Q* 
(Q'- -f 3) 3 



((e- 1)Q*^ -He-3) 



arctan Q* 



((e-l)Q*--H3e-)-l) 



arctan Q* 



(2Q'2 + 3)-3(n*^-t-l) 



arctan Q* 



Q.* 



